Hexatic, Wigner Crystal, and Superfluid Phases of Dipolar Bosons 



o 
o 



(N 



5h 



i 

C 

o 
o 



> 
in 
in 

o 

On 
O 



X 



Kaushik Mitra, C. J. Williams and C. A. R. Sa de Melo 
Joint Quantum Institute, University of Maryland, College Park, Maryland 20742, 
and National Institute of Standards and Technology, Gaithersburg, Maryland 20899 

(Dated: March 26, 2009) 

The finite temperature phase diagram of two-dimensional dipolar bosons versus dipolar interaction 
strength is discussed. We identify the stable phases as dipolar superfluid (DSF), dipolar Wigner 
crystal (DWC), dipolar hexatic fluid (DHF), and dipolar normal fluid (DNF). We also show that 
other interesting phases like dipolar supersolid (DSS) and dipolar hexatic superfluid (DHSF) are 
at least metastable, and can potentially be reached by thermal quenching. In particular, for large 
densities or strong dipolar interactions, we find that the DWC exists at low temperatures, but 
melts into a DHF at higher temperatures, where translational crystaline order is destroyed but 
orientational order is preserved. Upon further increase in temperature the DHF phase melts into 
the DNF, where both orientational and translational lattice order are absent. Lastly, we discuss the 
static structure factor for some of the stable phases and show that they can be identified via optical 
Bragg scattering measurements. 
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Arguably, one of the next frontiers in ultracold atomic 
and molecular physics is the study of ultracold het- 
eronuclear molecules such as KRb [ll, S Q , RbCs 
and NaCs which can be produced using Feshbach 
resonances observed in mixtures of two types of alkali 
atoms [1, 0,0- Thus, ultracold heteronuclear molecules 
consisting of Bose-Bose, Bose-Fermi, or Fermi-Fermi 
atom pairs offer many new op por tunities because of their 
internal degrees of freedom @,M EH, E 

[T^ such as their 

permanent electric dipole moment. 

The dipolar interaction between heteronuclear 
molecules is highly anisotropic in three dimensions hav- 
ing attractive and repulsive contributions. This makes 
it generally quite difficult to identify stable phases 
with good accuracy, fn the particular case of bosonic 
heteronuclear molecules (Bose-Bose or Fermi-Fermi) the 
attractive part of the dipolar interaction may lead to 
undesired instabilities of the dipolar gas. However, the 
situation in two dimensions (2D) can be quite different, 
and arguably more interesting, as the dipolar interaction 
can be made to be purely repulsive by the application of 
suitable static electric or microwave fields. In the case of 
bosonic dipolar molecules, several stable and metastable 
many body phases of 2D systems may be found. 

In this manuscript, we obtain the finite and zero- 
temperature phase diagram of bosons interacting via 
short-range repulsive interactions U and long-ranged 
dipolar interactions Ed in two-dimensions. For weakly 
repulsive values of U, and small values of Ed we find 
a dipolar superfluid phase of the Bcrczinskii-Kostcrlitz- 
Thoulcss (BKT) type, which upon increasing values of 
Ed becomes a dipolar Wigner crystal (DWC) forming 
a triangular lattice. Numerical evidence for the exis- 
tence of the DWC phase was obtained recently in Quan- 
tum Monte Carlo simulations of dipolar boson systems 
at zero temperature [Til [To]]. Here, however, we develop 
an analytical variational theory that accounts not only 



for the superfluid to DWC phase transition at zero tem- 
perature, but also for the finite temperature melting of 
the DWC into a dipolar hexatic fluid (DHF), where crys- 
talline translational order is destroyed but hexagonal ori- 
entational order is preserved. Further temperature in- 
crease leads to the melting of the hexatic phase into a 
dipolar normal fluid (DNF). We also find that the dipo- 
lar supersolid phase (DSS), exhibiting both superfluid 
and crystalline order, has lower (higher) energy than the 
DSF (DWC) as density or dipole energy increases, but is 
at least metastable, thus being accessible using thermal 
quenching. The DSS phase can also melt into a dipo- 
lar hexatic superfluid (DHSF). Lastly, we indicate that 
measurements of optical Bragg scattering can identify the 
DWC and DHF phases. 

Hamiltonian: To describe the phases of interacting 
dipolar bosons in 2D, we start with the continuum Hamil- 
tonian 
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where Vi oc = ^J2(i j) U5{ri —Tj), and the sum over 
indicate the sum over all pairs of molecules. The first 
term of H corresponds to the kinetic energy, the second 
to dipolar interactions, and the third (Vi oc ) to the local 
(short-range) interaction. 

We begin our discussion of different ground states by 
analyzing first the dipolar Wigner crystal (DWC) phase. 
In this phase the dipolar interactions are dominant in 
comparison to the kinetic energy and local energy terms, 
such that the system crystalizes into a triangular lattice 
in two dimensions (2D). Thus, our variational wavefunc- 
tion can be chosen to be 



*wc(ri,r 2 ,...) = -7=Y.Ii G ^ r ^ 

p(a) i 
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where {a^} forms a triangular lattice in 2D, and A is 
the normalization constant. The J2»(a ) * s a sum over 
all possible permutations P of {a} and the function 



G ai , a (ri) = e 



-(r<-a<) /2a' 



1 1 \pna is a normalized Gaussian 



centered at the lattice site a^. We define the separation 
between neighboring lattice sites to be a, and express 
the boson density as p = 2/(v / 3a 2 )- In addition, we in- 
troduce the dimensionless parameters rp = 2mDp 1 ^ 2 /h 2 
{rjj = 2mU /fi 2 ) as the ratio of the characteristic dipolar 
energy Ed = Dp z / 2 (local energy Ejj — Up) and kinetic 
energy K = h 2 p/2m. 

The variational wavcfunction described in Eq. J2]) is 
expected to be good for 3> 1, where the dipole-dipolc 
interaction is much larger than the kinetic energy, and 
can be used to compute analytically the average kinetic, 
local and dipolar energies in terms of the variational 
parameter a = a / 'a, where a is the Gaussian width and 
a is the lattice spacing. The dipolar Wigner crystal is 
only expected for < a < 1, where its average kinetic 
energy is (*|!f |*)dwc = NV3h 2 p/(4ma 2 ). The average 
local energy is (^\V loc \^) dwc = NSV^U pP(a)/ (2na 2 ) 
where P(a) = (e" 1 / 2 " 2 + e - 3 / 2 « 2 + e" 2 /" 2 ) and 

the average dipolar energy is (^Vdipl^dwc = 
N3V3KDp [P(a) + F(a)} /(2a) where F(a) 
(l + 9a 2 /2 + 225a 4 /8 + 3675a 6 /16). In general, 

P{a) and F{a) are represented by an infinite series in 
the variational parameter a, but the series is rapidly 
convergent for < a < 1, such that the first few terms 
are sufficient for the discussion of the total average 



The minimization of E, 



dwe 



energy E dwc = 
with respect to a leads to the minima illustrated in 
Fig. [TJ for ru = and a few values of rjj. For small 
values of rjj the transition is shifted towards lower values 
of td (not shown in Fig. [T]). 




1/a 

FIG. 1: (Color Online) Plots of the dipolar Wigner crystal 
energy -Edwc in units of the dipolar energy Ed versus the in- 
verse of the variational parameter a for ro = 15 (blue dashed 
line), td = 26 (red solid line), ro = 47 (green dotted line). 
The horizontal (black dot-dashed line) is the energy _Bd a f for 
a uniform superfluid. 

Superfiuid Phase: The superfluid phase is more easily 



described by writing our original Hamiltonian in Eq. ([T]) 
in second quantized notation, H = K + V, with the 
kinetic energy K = — j— / dr$ T (r)V 2< I>(r) and poten- 
tial energy V = \J draV$t( r )$t (r')V(r, r')$(r)$(r'), 
where V(t,t') = D/\r-r'\ 3 + U5(r-r'), and $ f (r) is the 
bosonic field operator which creates a dipolar molecule at 
position r. Describing the superfluid phase by the aver- 
age (4>t(r)) = ($(r)) = ^fp leads to the ground state 
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, by assuming that the 
3 is unscreened for length 



energy E dsf = N 

dipolar potential Vdi P (r) = D / 
scales r > ro, and screened to Vdi P (r) = D/tq for length 
scales r < r , where r « 0.9a. In Fig. [TJ the energy of 
the DSF phase is shown and compared with that of the 
DWC phase. 

Next, we discuss the melting of the Wigner crystal 
phase, which occurs in two stages. First the dipolar 
Wigner crystal melts into a hcxatic fluid, which docs 
not have translational order, but preserves rotational 
order. The melting occurs via the Kosterlitz-Thouless- 
Nelson-Halperin- Young (KTNHY) dislocation prolifera- 
tion mechanism. Second, the dipolar hexatic fluid trans- 
forms itself into the dipolar normal fluid by losing its 
rotational order at a higher temperature. 

Melting from Wigner Crystal to Hexatic Phase: To 
study the melting of the DWC phase, we need to calcu- 
late its elastic energy. This calculation can be performed 
by using a semiclassical approximation to the quantum- 
mechanical method for the calculation of elastic ener- 
gies [lij]. The essential idea is to stretch the many-body 
wavefunction \l/(ri,r2, ...) on each particle coordinate as 
Ti.a — * fi ta +e a pri/3 (where repeated indices indicate sum- 
mation), and expand up to second order in the strain 
tensor e a p. The resulting elastic energy is 



E* = \ J dv[2pe 2 ap + \e 2 aa \ , 



(3) 



where p = \b\/?>D / Aa 5 and A = 3p are the unrenormal- 
ized Lame coefficients, and the symmetric (rotation-free) 



strain tensor is e a /3( r ) 



du a (r) 



with u a (r) 



being the displacement from equilibrium position. 

We follow Ref. [l7j], and decompose the strain tensor 
e a p(r) into a regular (smoothly varying) cj) a p(r) and a 
singular (dislocation) d a p(r) contribution. Using k B = 1, 
the Hamiltonian for the smooth part can be written as 



T 



\ J V 2 ^ 



2 1 

aa I i 
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where p = pa 2 /T and A = Xa 2 /T. Additionally the 
Hamiltonian for the singular part can be written as 



T 



where b a is the alpha-component of the Burger's vector 
b defined by the contour integral of the displacement 



3 



field u around the dislocation: J du = ab(r). Also, the 
coefficient K = Kqo 2 /T, is related to the unrenormalizcd 
Lame constants through K = 4fi(fi + A)/(2/i + A), and 
E c = (i] + l)KT/8ir is the core energy associated with 
a dislocation of core diameter rya, while the interaction 
coefficient is 



a/3 



= Sa/3 hi 



r — r 



r — r 



(6) 



To obtain the critical temperature for the melting T m 
of the dipolar Wigner crystal we solve the renormaliza- 
tion flow equations for the elastic parameters 



djJL~ l 
dl 

dl 



3V(Z)e 



= 3V(0e K(0/8 "(F o (0+^i(0): 



of the Hamiltonian H lcg , where F n (l) = I n (K (I) / 8ir) , 
with I n (x) being the modified Bcsscl function of order n. 
These equations need to be solved sclf-consistcntly with 
the flow equations for the parameters 



dy(l) 
dl 

dK- 1 



2-^]y(l) + 2ny\l)e K ^F (l) 



of the dislocation Hamiltonian -ffdis- All flow equations 
are accurate to order y 3 , where y = e~ E "l T plays the role 
of the fugacity. The critical temperature Th is reached 
when 



K{T h ) 



2fi(T h ) + A(T fc ) 



16tt 



(7) 



Solving the RG flow equations, we get the melting tem- 
perature from the Wigner crystal to the hexatic phase to 
be Th = 0.05-Ed in the classical limit of m — ► oo. The 
melting temperature Th separating the DWC and DHF 
phases is shown in Fig. [2] 

Transition from Hexatic to Normal Phase: As the 
dipolar Wigner crystal melts at T/j, translational order 
disappears but oricntational order is preserved, with the 
emergence of the hexatic order parameter ^q{t) = e 6l9 ( r \ 
The bond-angle field 0(r) between the location of the 
center of masses of the dipolar bosons (heteronuclcar 
molecules) is related to the displacement field u(r) by 
6(r) = [d x u.y(r) — d y u x (r)] /2. The elastic energy for 
such a situation can be also obtained using a semiclas- 
sical approximation to the method of Ref. [IB], where 
each particle coordinate in the many-body wavefunction 
^(ri, T2, ■•■) is locally rotated ri <a — > M a pri,p, where 
M a /s is a local rotation matrix (tensor) . For the triangu- 
lar lattice considered here the elastic energy becomes 
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FIG. 2: Finite temperature phase diagram of T/Ed versus ro 
showing the dipolar superfluid (DSS), dipolar Wigner crys- 
tal (DWC), dipolar hexatic (DHF) and dipolar normal fluid 
(DNF) phases, phases. 



where Tq w 2E c is the phase stiffness of the hexatic 
phase, which is directly related to the dislocation core 
energy E c = l.lKT/8w. The critical temperature for 
the disappearance of hexatic order and the emergence of 
the normal phase T„ is then determined by the RG flow 
of the 2D-XY model [lq . which leads to the condition 
^6{T n ) = 72T n /ir. The algebraic decay of orientational 
order is then destroyed at temperature T n by the dis- 
sociation of pairs of ±7r/3 disclinations, which play the 
role of vortices and anti- vortices of the standard 2D-XY 
model. In Fig. [21 the critical temperature T n separating 
the phases DHF and DNF is shown as a function of re. 
When m — ► oo, we obtain T n — > Q.WEd which is larger 
than the melting temperature of the dipolar Wigner crys- 
tal. 
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FIG. 3: Phase diagram showing possible dipolar supersolid 
(DSS) and dipolar hexatic superfluid (DHSF) phases. 

Supersolid Phase: Starting from the Hamiltonian de- 
fined in the superfluid section we seek a variational solu- 
tion 

($(r)) ss = $ ss (r) = N — • x — (9) 

i 

for the supersolid phase, where the Gaussian functions 
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G ai ,<T ss (r) form a triangular lattice of side a ss . The 
normalization condition is J dr\Q ss (r)\ 2 /V = p, while 
(T S s = P&ss is the gaussian width and /3 is the correspond- 
ing variational parameter. The supersolid order param- 
eter $ ss (r) describes a non-uniform superfluid with both 
off-diagonal long-range order due to broken U(l) sym- 
metry and diagonal long-range order due to the discrete 
lattice symmetry. When p s f 3> p ss , the energy is es- 
sentially that of a superfluid, and when p s f <C p ss , the 
average kinetic energy is K dss = NV3h 2 p/(4m/3 2 ), while 
the total potential energy is Vdss = V\ + V2, where Vi = 
N3V3UpP(f3)/(27rP 2 )+NV3Up/8irP 2 is the local poten- 
tial energy and V 2 = NZVSKDp [P(J3) + F(J3)] /(2a ss ) + 
NiJ?,Dp/Aal s is the dipolar energy. 

We find that the total energy Ea ss = -Kdss + Kiss of the 
dipolar supersolid (DSS) phase can be lower than that 
of the superfluid phase due to the contributions coming 
from the dipolar interaction, but E^ ss is always higher 
than the total energy i?dwc of the dipolar Wigner crystal 
for the same values of parameters. Thus, the transition 
from the superfluid to the supersolid phase does not oc- 
cur within our variational ansatz, even if we include ad- 
ditional correlations via Jastrow factors [la [. However, 
this transition can not be excluded, because only a few 
classes of variational wavefunctions have been explored. 
Additionally, we find that the supersolid phase is at least 
metastable, since E ss has a minimum that can be reached 
upon thermal quenching, and possibly also melted into a 
metastable dipolar hexatic superfluid (DHSF). A phase 
diagram showing possible DSS and DHSF phases is pre- 
sented in Fig. [3J 

Experimental characterization of the various phases: 
The various phases proposed here can be characterized 
by the measurement of their density-density correlations, 
as reflected in the structure factor S'(q) = (n(q)n(— q)), 
where n(q) is the Fourier transform of the density opera- 
tor n(r), and the symbol (...) corresponds to both quan- 
tum mechanical and thermal averaging. 



a) • 



FIG. 4: Bragg scattering patterns near the first reciprocal 
vectors for a) the dipolar Wigner crystal phase and for b) the 
dipolar hexatic fluid phase. 

The most dramatic effects in the structure factor are 
found in dipolar Wigner crystal phase where S(q) ~ 
|q — G|~ 2+ly ( G,T ) revealing the power law behavior char- 
acteristic of two dimensions in the vicinity of the recip- 
rocal lattice (Bragg) vectors G. The first Bragg vector 



has magnitude |Gi| = 47r/\/3a, and the Bragg scattering 
pattern then exhibits 6-fold symmetry below the melt- 
ing temperature T m as shown in Fig. [4^. The exponent 
v(G,T) = (|Gi|/47r)(T//i„)(3M. + A»)/(2/i* + A*) is re- 
lated to the renormalized Lame coefficients /1* and A*. 
This reflects the power law decay of the correlation func- 
tion C(G,R) — |R|" l/(G ' T) , which is the Fourier trans- 
form of 5(q). The profile of S(q) for the dipolar hexatic 
phase corresponding to a melted dipolar Wigner crystal 
with orientational order is characterized by the hexago- 
nal pattern shown in Fig. HJd. 

Conclusions: We have described some of the possible 
phases of dipolar bosons in two dimensions, including 
superfluid, supersolid, hexatic superfluid, Wigner crys- 
tal, hexatic and normal fluids. Within our variational 
approach we have concluded that the supersolid and hex- 
atic superfluid phases may have lower free energy than 
the superfluid phase in some region of the phase dia- 
gram, but do not seem to have lower free energy than the 
Wigner crystal or hexatic fluid in the same region. How- 
ever, the supersolid and hexatic superfluid phases are at 
least metastable, and thus may be reached via thermal 
quenching, and probed via Bragg scattering. Further- 
more, we have shown that dipolar Wigner crystal does 
not melt directly into a normal fluid, but presents a two- 
stage melting first into a hexatic phase which preserves 
orientational order, and then into a normal fluid upon 
further increase in temperature. Finally, we indicated 
the experimental signatures of stable phases with trans- 
lational or orientational order in a Bragg spectroscopy 
measurement that detect the static structure factor. 
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